Transport regimes in surface disordered graphene sheets 
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We investigate the size scaling of the conductance of surface disordered graphene sheets of width 
W and length L. Metallic leads are attached to the sample ends across its width. At the Dirac 
point, E = 0, the conductance scales with the system size as follows: i) For constant W/L, it remains 
constant as size is increased, at a value which depends almost linearly on that ratio ; this scaling 
allows the definition of a conductivity value that results similar to the experimental one. ii) For 
fixed width, the conductance decreases exponentially with length L, both for ordered and disordered 
samples. Disorder reduces the exponential decay, leading to a higher conductance, iii) For constant 
length, conductance increases linearly with width W , a result that is exclusively due to the tails of 
the states of the metallic wide contact, iv) The average conductance does not show an appreciable 
dependence on magnetic field until fields such that the flue per unit cell approaches the quantum 
unit. Away from E = 0, the conductance shows the behavior expected in two-dimensional systems 
with surface disorder, i.e., ballistic transport. 

PACS numbers: 73.63.Fg, 71. 15. Mb 
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Introduction. The electronic transport in atomically thin 
graphene samples is a subject of great current interest 
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The scaling with the sample dimen- 
sions |4| suggest a diffusive behavior, with a universal 
conductivity at the lowest carrier concentrations [1, @. 
The limit of low concentrations is difficult to analyze 
theoretically, as the Fermi wavelength becomes compa- 
rable to the separation between scatterers, and even to 
the sample size. An analysis based on the Born approxi- 
mation 8| leads to a universal conductivity at low tem- 
peratures, although its value is somewhat smaller than 
the one observed experimentally. The approximations 
involved in this approach, however, are expected to fail 
at the lowest concentrations. Field theoretical arguments 
Ei U | suggest the existence of a localized regime in 
the limit of zero temperature and zero carrier concentra- 
tion. At zero doping, clean graphene systems show an 
unusual scaling of the conductance on sample size, con- 
sistent with diffusive behavior [l2| ■ This pseudo-diffusive 
behavior has also been found in SNS junctions [l3| and 
graphene bilayers 14\. 



In this work, we numerically study the electronic trans- 
port in surface disordered graphene sheets both at finite 
dopings and in the limit of zero carrier concentration. As 
bulk disorder in graphene sheets seems to be rather low, 
we focus on the effects of rough edges, with disorder con- 
centrated at the surface of the system. Our results show 
that the pseudo-diffusive regime identified in [l^] persists 
in the presence of disorder, namely, near the band center 
the conductance is proportional to the sheet width and 
almost inversely proportional to its length. Although our 



results fit apparently the requirements of diffusive scaling 
in 2D, a closer look reveals important differences. Cer- 
tainly, when plotting the conductance of samples of fixed 
width as a function of the sample length, an exponential 
decrease is obtained, that is, the standard result for a 
quasi ID system with any kind of disorder. The remark- 
able thing in graphene is that this exponential decrease 
survives in ordered samples. Actually, the presence of 
disorder slows down the exponential decay. These results 
can be rationalized in terms of transmission mediated by 
evanescent waves generated at the metallic leads. On 
the other hand, the conductance for fixed length is pro- 
portional to the sample width at all energies, a behavior 
that does not distinguish between diffusive and ballistic 
regimes. Finally, for sufficiently high carrier concentra- 
tions, the conductance shows the ballistic behavior ex- 
pected in 2D systems with surface disorder, namely, in- 
crease with the system size for constant W/L, linear in- 
crease with W for constant L, and exponential decrease 
with L for constant W. 

Methods: graphene samples and conductance calcula- 
tions. We describe the valence and conduction bands 
of graphene by a tight binding hamiltonian with nearest 
neighbor hoppings only: 



4 C 3 



(1) 



where sites i and j denote the nearest neighbor nodes in 
the honeycomb lattice. The low energy electronic spec- 
trum, |ej-| Ct, can be approximated by the Dirac equa- 
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FIG. 1: (Color online) Scaling of the conductance of stripes 
of different widths, W, as function of length, L. 



tion: 



±v F \k\ 



(2) 



where vf = (3td)/2, and d is the distance between sites 
in the honeycomb lattice. 

Samples with surface disorder were produced by ran- 
domly removing sites at the sheet edges. The leads were 
simulated by a purely imaginary selfenergy independent 
of energy, that was attached across the sample width. 
In the interpretation of the numerical results for surface 
disordered sheets, perfect samples without and with An- 
derson disorder at its edges lH, [lj| were also investi- 
gated . The latter was introduced by randomly sorting 
the orbital energies at the surface sites within the range 
[-A,A]. 

The conductance was calculated by means of an ef- 
ficient implementation of Kubo's formalism 17]. For a 
current propagating in the x-direction, the static electri- 
cal conductivity is given by: 



-2|-]Tr 



(M x )Im G(E) (hv x )Im G(E) 



(3) 



where the velocity (current) operator v x is related to 
the position operator x through the equation of motion 



ttv 



H,x 



H being the Hamiltonian. G(E) is the 
Green function of the system with the leads already in- 
corporated. All results include the spin degeneracy. 
Results. Fig.pQ shows typical results for stripes without 
and with disorder. At E = 0, the scaling of the con- 
ductance in clean samples as G — 4e 2 / (tin) x VF/L[l3] is 
already obtained with high accuracy in not too wide sam- 
ples. At higher energies, E = 0.5, the conductance be- 
comes ballistic. This pseudodiffusive regime, was already 
analyzed in clean systems in [12j . In a clean square sys- 
tem, the incoming channels can be characterized by the 
transverse momentum, k y . The electronic spectrum of a 
graphene stripe at finite transverse momentum shows a 
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FIG. 2: (Color online) Top panel: Scaling with length L of 
the conductance in samples with width over length ratios and 
surface disorder. Bottom panel. Scaling with length of stripes 
of fixed width. The points labelled metallic correspond to a 
stripe with a subband crossing at E = 0. All other samples 
have finite size gaps near E = 0. 



gap for — i>F|fc«| < e < wp|/cj,|. Hence, states with trans- 
verse momentum k y decay away from the boundaries as 
e~\ ky \ x , and lead to a transmission Tk v oc e~ 2 \ k y\ L , where 
L is the length of the system. The number of channels 
scale as the width of the system, W. In sufficiently large 
systems, the sum over channels can be replaced by an in- 
tegral over k v leading to a conductance G which scales as 
G oc WL^ 1 . Away from E = the conductance increases 
linearly with the system size. This is the expected bal- 
listic behavior of a quantum billiard with either surface 
disorder or with an amount of defects proportional to L 

urn. 



(a defect concentration decreasing as 1/L) 

The scaling of the conductance with stripe length in 
systems with surface disorder and different W/L ratios 
is shown in the upper panel of Fig. [2j. The conduc- 
tance at E = 0.1 deviates very slightly from ballistic 
behavior. Disordered samples show a length indepen- 
dent conductivity close to that estimated analytically in 
the clean linmit 12| , and this regime is attained even for 
widths smaller than the length. The conductance scales 
exponentially with L. In a clean systems, this behavior 
arises from the existence of minigaps separating subbands 
with well defined periodicity in the transverse direction. 
The corresponding decay length shorter in clean samples 
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FIG. 3: (Color online) Conductance (in units of the conduc- 
tance quantum) through surface disordered graphene samples 
of size W x 24d versus the sample width W (in units of the 
C-C distance a). Results for three energies are shown. 



than in disordered samples, leading to an enhancement 
of the conductance in semiconducting disordered samples 
at E = 0. 

The enhancement of the conductance in the localized 
regime in disordered samples is probably related to the 
formation of resonances at E = near defects. The ex- 
istence of these resonances has been well established at 
edges [l9l [2fj| . cracks [2lj|. and vacancies [22j|. After the 
first version of this paper was posted, related resonances 
where discussed in the continuum limit [23| . Note that, 
in addition to E = resonances induced by disorder, 
long wavelength modulations of the chemical potential 
will move one of the edges of the gap at fixed parallel 
momentum towards E — 0, reducing the decay length 
and enhancing the conductance. 

The scaling of the conductance with the width of a 
sheet of constant length is shown in Fig. [3] The numeri- 
cal results clearly indicate that the conductance increases 
linearly with the sample width with a slope that depends 
on the concentration of carriers (or the energy) . This re- 
sult is characteristic of both ballistic and diffusive behav- 
iors in 2D, and cannot therefore be used to discriminate 
the transport regime in this case. A remarkable feature 
of the results shown in Fig. [3] is its very low dispersion. 
This could be understood by noting that the increase in 
conductance is exclusively due to an increasing number 
of metallic tails through the bulk of the graphene sheet, 
and therefore weakly sensitive to surface disorder. 

We show in Fig. [4] the dependence of the conductance 
on magnetic field for different disorder realizations [2^ |. 
The area of the sample, 40d x 40d, and 62 unit cells, is 
such that one flux unit through it is equivalent to ap- 
proximately 660 Teslas. Hence, the magnetoresistance 
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FIG. 4: (Color online) Conductance (in units of the conduc- 
tance quantum) as function of the total magnetic flux through 
the graphene lattice in 40d x 40d clusters and E = 0.02, with 
disorder at the edges. A magnetic field of 1 Tesla corresponds, 
approximately, to 0.04 flux units through the cluster. 



for fields attainable in the laboratory is negligible. This 
result is consistent with semianalytical calculations using 
a continuum model in the clean limit [2^]. The magneto- 
resistance shows oscillations when the flux per unit cell 
is of order 0.05. At very high fields, the conductance 
becomes of order of one quantum unit, and it shows no 
dependence on disorder realization. Note that, in the 
regime studied here, the quasiclassical arguments used 
when discussing either weak localization or weak antilo- 
calization effects in graphene cannot be used SEES. 

The study of the conductance distribution for samples 
of an approximate square shape at an energy E — O.OOli 
very close to the Dirac point further explains the role 
played by the metallic tails in the conductance behav- 
ior. Fig. [5] shows that conductance is larger than 
w 0.37(2e 2 //i) for this geometry (the tails contribution 
that is minimally affected by surface disorder) and fluc- 
tuates below 1 as it does in a standard quantum billiard 
in the case of point contacts. The existence of an abrupt 
upper cutoff resembles the case studied in[29(. Neverthe- 
less, a significant difference is clear; while the upper con- 
ductance limit is due to the incidence of only one channel 
in the billiard case, it is due to the intrinsic small number 
of channels (0 or 1) of graphene near the band-center. 
Concluding Remarks. The numerical calculations of the 
conductance through surface disordered graphene sheets 
presented in this work reproduce the quasi-diffusive be- 
havior foun d by other authors in ordered graphene at the 
Dirac point 12, 13, 14, 3(J. Specifically, the conductance 
remains constant when the size of the system is increased, 
as opposed to the linear increase with the system size 
found at any other energy. However, we found this be- 
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FIG. 5: Conductance distribution obtained at E — O.OOlt for 
ia set of 10,000 graphene randomly generated samples of size 
24d x 24d. Surface disorder is restricted to a Ad fringe around 
the sample surface. Good metallic contacts are attached at 
opposite sides of the sample. 

havior only when the size of the system is increased in 
such a way that the width to length ratio is kept constant. 
A value of w 0.75(2e 2 //i) is obtained for this pseudo- 
conductivity, not that far from the experimental one. 

The conductance of the stripes calculated here changes 
qualitatively as the aspect ratio of the sample is var- 
ied. In stripes where the length is much larger that the 
width, the pseudodiffusivc behavior described above is 
replaced by an exponential decay with length, a sign of 
localization. In this regime, disorder leads to longer de- 
cay lengths, probably due to the formation of resonances 
at the Dirac point, E = 0. 

Finally, we have also shown that a magnetic held with 
a magnetic length much larger than the lattice spacing 
does not change appreciably the conductance near the 
Dirac point. 
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